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Q . Abstract 
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1-Q \ Within an effective field theory framework, we obtain an expression, with 

X.' 0(1/?TT-^) accuracy, for the energies of the gluonic excitations between heavy 

(— I \ quarks, which holds beyond perturbation theory. For the singlet heavy quark- 

antiquark energy, in particular, we also obtain an expression in terms of Wilson 

kS ' loops. This provides, twenty years after the seminal work of Eichten and 

Feinberg, the first complete expression for the heavy quarkonium potential 
up to 0{\/m?) for pure gluodynamics. Several errors present in the previous 
literature (also in the work of Eichten and Feinberg) have been corrected. We 
also briefly discuss the power counting of NRQCD in the non-perturbative 
regime. 
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I. INTRODUCTION 

The measured spectroscopy suggests that the charm and bottom quark masses are large 
enough to consider their heavy-quark-antiquark bound-state systems (generically denoted 
as heavy quarkonia: ip, T, Be-, ...) as non-relativistic (NR). These systems are, therefore, 
characterized by, at least, three widely separated scales: hard (the mass m of the heavy 
quarks), soft (the relative momentum of the heavy-quark-antiquark |p| ~ mv, v <^ 1), and 
ultrasoft (the typical kinetic energy E ~ mv"^ of the heavy quark in the bound-state system). 
Inspired by this NR behaviour, the investigation of heavy quarkonia has been traditionally 
performed by all sorts of potential models, where an ansatz potential is introduced in a 
Schrodinger equation (for some reviews see 0^)- The phenomenological success of these 
suggests that, to some extent, a potential picture may, in fact, be appropriate and justified 
from QCD. This triggered the attempts to derive these potentials from QCD by relating 
them to Wilson loops. These standard derivations used an expansion in 1/m (also named 
adiabatic or Born-Oppenheimer approximation). However, a full derivation of the potential 
from QCD, as well as a study of the validity of the potential picture itself, was not done so 
far in the non-perturbative regime, where most of the heavy quarkonium spectrum lies. It is 
the aim of this paper to explicitly derive the complete non-perturbative 1/m? QCD potential 
for pure gluodynamics within an effective field theory framework Q,^, where higher order 
potentials in 1/m and non-potential effects could also be incorporated in a systematic way. 

Since the derivation of the potential has a long story, it may be useful to summarize 
its main steps. The expression for the leading spin-independent potential, of 0{l/m^), 
corresponds to the static Wilson loop and was derived and discussed in the seminal works 
of Wilson and Susskind @J^. Expressions for the leading spin-dependent potentials in the 
1/m expansion, of 0{l/m?), were given in Refs. [PHTOH. The procedure followed in these 
works proved to be very difficult to extend beyond these leading-order potentials. Indeed, 
the first attempts [ll| , using tools similar to those in Ref . [Q , failed to obtain suitable finite 



expressions. In Ref. []T^, a new method to calculate the potentials was proposed, where 
new spin-independent (some of them momentum-dependent) potentials at 0{l/m?') were 
obtained. In these original works, the obtained potentials did not correctly reproduce the 
ultraviolet behaviour expected from perturbative QCD (the hard logs ~ \ogm). This was 
first implemented in the framework of QCD effective field theories, for both spin-dependent 
and spin-independent potentials, in [^^-|T5|. At that point, the obtained set of potentials at 
0{l/m?) seemed to be complete and the timely study of the different Wilson loop operators 
describing the non-perturbative dynamics of the potentials started. For instance, a lattice 
study was performed in [Q and a study in the framework of QCD vacuum models was done 



m 



Nevertheless, this view has been recently challenged in Ref. where: i) a systematic 
study of the potentials has been started within an effective field theory framework: potential- 
NRQCD (pNRQCD) 0], and ii) the 0(l/m) potential, previously missed in the literature, 
has been calculated. It is the aim of this paper to explain in more detail the Hamiltonian 
formalism, sketched in Ref. 0, and to compute the 0{l/m?) potentials. The formalism 
appears to be quite powerful and suitable to obtain the quarkonium potentials and the 
energies of any gluonic excitation at any finite order in 1/m. A similar idea, but in the 
Coulomb gauge and only for the leading spin-dependent quarkonium potentials, has also 



been used in |T^. We will give an expression in terms of quantum-mechanical corrections 
to the energies of the gluonic excitations between static quarks, valid for all the gluonic 
excitations up to 0(l/m^). For the quarkonium state (the ground state), we will express 
our complete 1/m? result in terms of Wilson loops eventually calculable on the lattice or by 
means of QCD vacuum models, concluding in this way an ideal journey started over twenty 
years ago. 

The theoretical framework of our work is NRQCD [|l^ and pNRQCD, suitable effective 
field theories for systems made up by two heavy quarks. NRQCD has proved to be extremely 
successful in studying heavy quark-ant iquark systems near threshold. It is obtained from 
QCD by integrating out the hard scale m. It is characterized by an ultraviolet cut-off much 
smaller than the mass m and much larger than any other scale, in particular much larger 
than Aqcd- This means that the matching from QCD to NRQCD can always be done 
perturbatively, as well as within an expansion in 1/m ]19|,20|. The Lagrangian of NRQCD 



can also be organized in powers of 1/m, thus making explicit the non-relativistic nature of 
the physical systems. So far, NRQCD and pNRQCD have only been studied in detail in the 
perturbative situation pT|J^. 



By integrating out degrees of freedom with energies larger than mf^, one is left to a 
new effective field theory called pNRQCD where the soft and ultrasoft scales have been 
disentangled and where the connection between NRQCD and a NR quantum-mechanical 
description of the system can be formalized in a systematic way. pNRQCD has two ultraviolet 
cut-offs, Ai and A2. The former fulfils the relation mv"^ <^ Ai <^ mv and is the cut-off of the 
energy of the quarks, and of the energy and the momentum of the gluons, whereas the latter 
fulfils mv -C A2 <^ m and is the cut-off of the relative momentum of the quark-antiquark 
system, p. In the non-perturbative situation (we understand by non-perturbative a typical 
situation where mv ~ Aqcd, i-e. where the potential cannot be computed perturbatively), 
we will assume that the matching between NRQCD and pNRQCD can be performed, as 
in the perturbative case, order by order in the 1/m expansion. We will present, for the 
general situation Aqcd < '<^v, the matching of NRQCD to pNRQCD at 0(l/m^) for the 
singlet sector (to be defined later). This will prove to be equivalent to computing the 
heavy quarkonium potential that we can now derive from QCD by a systematic procedure. 
Moreover, the expression for the potential that we obtain will also be correct at any power 
in Os in the perturbative regime. 

A pure potential picture emerges in pure gluodynamics under the condition that all 
gluonic excitations have a gap larger than mf ^. Extra ultrasoft degrees of freedom such 
as hybrids and pions can be systematically included and may eventually affect the leading 
potential picture (as ultrasoft gluons in the perturbative regime 0). 

In this paper we consider the general situation of particles with different masses. There- 
fore, our results, besides to the traditional Q-Q systems, may be applied to the Be system. 



which, after its recent discovery by the CDF collaboration [22|, has received a lot of attention 



in theoretical investigations [23 



The paper is organized in the following way. In section ^ we introduce NRQCD up to 
0(l/m^). In section |T|, using a Hamiltonian formulation of NRQCD, we explicitly calculate 



up to 0{l/m?) the energies of the gluonic excitations between heavy quarks. In section |V] 
we define what pNRQCD will be in the present context. In section we write the heavy 
quarkonium potential up to 0(l/m^) in terms of Wilson loops and compare with previous 



results. In section |VT| we discuss the power counting of pNRQCD in the non-perturbative 
regime and in section [VI I| we give our conclusions and outlook. 



II. NRQCD 

After integrating out the hard scale m, one obtains NRQCD [|18[. Neglecting operators 



that involve light quark fields ||2J], the most general NRQCD Lagrangian (up to field re- 
definitions) for a quark of mass mi and an antiquark of mass m2 up to 0(l/m^) is given 
by: 




^_, [D-,E] , . (2) ^-[Dx.E] ] 

67712 01712 J 

mim2 77111712 mim2 mi7n2 

--G^,G + |^_ + —J G^,/; G + |^_ + —J (7A,eG^,G^„G,„, 

where il) is the Pauli spinor field that annihilates the fermion and x is the Pauli spinor field 
that creates the antifermion, iD^ = ido — gA°, iD = iV + gA, [D-, E] = D ■ E — E ■ D and 
[Dx, E] = D X E — E X D. This Lagrangian is sufficient to obtain the 0(l/?7i^) potentials. 



The coefficients cp, cd, cs, d2 and d^ can be found in Ref. |]T9l and dij {i,j = s,v) in |2^ for 



the MS scheme. 

Some words of caution are in order here. Even if the above matching coefficients have 

been computed using dimensional regularization and the MS scheme, there could still remain 

some ambiguity depending on the different prescriptions for the e^^'' tensors and the definition 

of the Pauli matrices cr. For instance, the use of a scheme where the e^^'' only takes values 

for dimension equal to three ('t Hooft-Veltmann-like scheme) in the computation of Ref. 

2 
SI| would change the value of dyy as dj^^ — > rf^^,, where D is the number of space-time 



dimensions. One should therefore be careful and make sure that the matching coefficients one 
is working with really are computed in the same scheme. A deep study of these ambiguities 
in the framework of NRQCD remains to be done. This may be specially important for higher 
order calculations. See also Refs. P3|,|2B[|, where the authors have to deal with equivalent 



problems. 

We are interested in the Hamiltonian of the above Lagrangian. The construction of the 
Hamiltonian of one effective (non-renormalizable) Lagrangian may be complicated (for a 
related discussion we refer to |^); in particular because there are higher time derivatives 
acting on the different fields. In order to get rid of those at 0(l/?7i^) we have to eliminate 
the term G't^D'^G'^^^ from the Lagrangian. This can be achieved by a field redefinition as 
follows. We consider the field redefinition of the gluon field (c ~ 1/m?): 

A^^A^ + c[D",G,,,]+0(c2), (2) 



where c is real. This transformation preserves the gauge transformation properties and the 
hermiticity of the A^ field. Eq. (^ produces the following change in the gluon Lagrangian 
(at the order of interest): 

- \g;^G-^'' -. -\gI,G-^^ - ^-G^D^C^^^ - cgU^Gl^Gl^Gl^ + 0{^) . (3) 

We can therefore cancel the GD^G term by fixing 

mf 1712 

This changes the value of (is to d'^: 

rf«' = 4') - 2rf« , 4^)' = df^ - 2df . (5) 

Let us now see the modifications that the above field redefinition will produce in other sectors 
of the theory, in particular, in the heavy fermion bilinear Lagrangian. Since we have the 
following change for the Do covariant derivative that appears at 0(l/m°) 

iD^^iD^-cg\D-,^\, (6) 

the matching coefficients change, at 0(l/'m^), as (all the others remain unchanged): 

c«' = c« - 164^) - 1644'^ , cg^' = cg^ - 164') - 16^4'^ • (7) 

In summary, eliminating the term G"" D'^G"'^^ ^ up to order \jvr? ^ is equivalent to the redef- 
inition of the matching coefficients ^3 — > dg and cd — * c^ found above. We will assume this 
field redefinition in the following. 

III. GLUONIC EXCITATIONS IN A HAMILTONIAN FORMULATION 

The Hamiltonian associated to the Lagrangian (|ip is, up to order 1/m^, 

mi 1712 mf 1712 mim2 

i7(°) = f d^x- (n'^n'^ + B°B'^) , (9) 

if (i'°) = -i / rf^t (d^ + gc?a . b) ^, i7(«'i) = i / rf^ (d^ + gcPa ■ b) x, (10) 

H^'''^ = Jd^^^P^i^-cP'g^^ - 4^) , " • ["g"' "^J ]^-Jd^^d^^^U.Gl^Gl^G:^, (11) 

^(0,2)^^(2,0)(^^^.1^2), (12) 

H^'''^ = -Jd'^ {dss^^^xh +ds,^^(T^x^(Tx +t/,,V^^T>x^T'^x +4.^^TV^x^TVx) , (13) 



where 11" is the canonical momentum conjugated to A'^ and the physical states are con- 
strained to satisfy the Gauss law: 

D ■ n'^lphys) = giiJ^T'^iJ + x^T«x)|phys)- (14) 

Since 11" = E" + 0(l/m^), in Eqs. ([TT| - |T2D and in the rest of the paper, we will use the 



chromoelectric field instead of the canonical momentum where, to the order we are interested 
in, it does not affect our results. 

A. The static limit 

We are interested in the one-quark-one-antiquark sector of the Fock space. In the static 
limit the one-quark-one-antiquark sector of the Fock space can be spanned by 

|n;xi,X2)(°)=^T(x,)^t(x2)|n;xi,X2)(°\ Vxi,X2, (15) 

where |n; xi, X2)^°-' is a gauge-invariant eigenstate (up to a phase) of H^^\ as a consequence 
of the Gauss law, with energy ii^^°^(xi,X2). For convenience, we use here the field Xc(x) = 
icr^X*(x)) instead of x(x)) because it is the one to which a particle interpretation can be easily 
given: it corresponds to a Pauli spinor that annihilates a fermion in the 3* representation of 
color SU{?>) with the standard, particle-like, spin structure. \n] xi, X2)^°-' encodes the gluonic 
content of the state, namely it is annihilated by Xc(x) and '?/'(x) (Vx). It transforms as a 
3x1 ® 3*2 under colour S'f/(3). The normalizations are taken as follows 

(°)(m;xi,X2|n;xi,X2)^°^ = (^nm, 

(°)(m; xi, X2|n; yi, y^)^'^ = ^n^^^^^^xi - yi)<5(='Hx2 - y2) • 

We have made it explicit that the positions xi and X2 of the quark and antiquark respectively 
are good quantum numbers for the static solution |n; Xi, X2)*-°'', whereas n generically denotes 
the remaining quantum numbers, which are classified by the irreducible representations 
of the symmetry group D^oh (substituting the parity generator by CP). We also choose 
the basis such that T|n; Xi,X2)'-°'' = |n; Xi,X2)*-°^ where T is the time- inversion operator. 
The ground-state energy E^ ■^(xi,X2) can be associated to the static potential of the heavy 
quarkonium under some circumstances (see Sec. [TVl ). The remaining energies ii^^°^(xi,X2), 
n 7^ 0, are usually associated to the potential used in order to describe heavy hybrids or heavy 
quarkonium (or other heavy hybrids) plus glueballs (see Sec. |VD. They can be computed on 
the lattice (see for instance ||2^). Translational invariance implies that £'^°''(xi, X2) = E^\r), 
where r = xi — X2. 



B. Beyond the static Umit 

Beyond the static limit, but still working order by order in 1/m, the normalized eigen- 
states, |n; xi, X2), and eigenvalues, -En(xi, X2; pi, P2), of the Hamiltonian H satisfy the equa- 
tions 



i7|n;xi,X2) =|d^x[d'x',\n■^[,^',)En{^[,^',■p[,p',)6^^\^[-^^)6^^\^',-^2), (16) 

(m; xi, X2|n; yi, y2) = 5nm5^^\xi - yi)5(^)(x2 - y2). (17) 

Note that the positions xi and X2 of the static solution still label the states even if the 
position operator does not commute with H beyond the static limit. We are interested in the 
eigenvalues En, which should be understood as operators (instead of numbers, even though 
we call them energies). This will match the operator interpretation within a quantum- 
mechanical formulation that we will give to them in pNRQCD in the next section. In 
particular, we will see that Eq corresponds to the quantum-mechanical Hamiltonian of the 
heavy quarkonium (in some specific situation). The other energies. En for n > 0, are related 
to the quantum-mechanical Hamiltonians of the heavy hybrids or heavy quarkonium (or 
other heavy hybrids) plus glueballs. 

Since the derivation of the corrections to En may not be familiar to the reader, since they 
are operators, we explain it in some detail. We will work in the same way as in standard 
quantum mechanics, but taking into account the fact that they are operators. Analogously 
to standard quantum mechanics, we define a state |n;xi,X2) such that 

H\n; xi, X2) = / (fx[d^x'2\n; x.[, X2)£'„(xi, X2; p[, P2)6''^\x.[ - Xi)5^^)(x2 - X2), 
(°)(n;xi,X2|n;yi,y2) = 5^'\^^ - y,)5^'\^2 - 72)- 
Splitting the Hamiltonian as H = Hq + Hj we have 

|n;xi,X2) = |n;xi,X2)(°) + ,^. — - J2 /c^Mrf^^lm; x'^, x'2)(°)(°)(m; x'^, x'sl 

En [X) — -O *■ -' nij^n 

xi^Hi\n; xi, X2) - J d^x'^d'^x'^ln; x[, x'2) A^„(x;, x'2; p'l, P2)5^''\A - ^i)5^'^\A - ^2)}, 

and 

AEn{xi,X2;pi,p2)5^^\xi-yi)5^^\x2-y2) = ^°^(n; xi,X2|i/7|n; yi, ys). 

From these formulas we can obtain En order by order in the expansion parameter of Hj. 
Moreover |n; Xi, X2) and En are given by 

|n;xi,X2) = y'rfVirfV2|n;x;,x'2)iV„-i/2(^/^^/.p/^p/)^(3)(^/ _ ^^)s^^)^^>^ _ ^2), 
and 

J-^n ^^n ^n-^^n ' 

where 

(n;xi,X2|n;yi,y2) = A^„(xi,X2; pi, P2)5(^)(xi - yi)5(^)(x2 - y2). 

By using the above results, we get for En up to 0{l/m'^): 



Eni^u X2; pi, P2)5(=')(xi - x;)5(=')(x2 - x'2) = i?i°)(xi, ^2)S'^'\^^ - ^[)6'^'\^2 - x'2 



2) 



t,,, ,i/(l'0) //(O'l) i/(2,0) ^(0,2) ^(1,1) 

+ (°) n; xi, X2 + + — + — + n; x'^, x^ ^^^ 

TJii m2 rrii 1712 11111712 



1 r f7{l>0) f7(0,l) 

- ^E /AiA2(°^(n;xi,X2|^i— + ^|k;yi,y2)(°) 



k^n 



x(°)(k;yi,y2 



^(1,0) ^(0,1) 

+ 



mi 



1712 



n;x;,X2) 



'\(0) 



X 



.(0), 



+ 



^r(yi,y2) - i5;r(xi,x'2) E^^'iy^^y^) - i5;r(xi,x2 



(0), 



^(0), 



The expansion of En in inverse powers of the mass can be organized up to 0(1/ 171"^) as 
follows: 



/?(i,o) jp{o,i) jp{2,o) s;(o,2) s;(i,i) 
mi m,2 m^i ??i2 m,im,2 



(19) 



From Eq. ([TsD and Eqs. (|ToD-(|T^), by using the identities (here and in the rest of the paper, 
if not explicitly stated, the dependence on Xi and X2 is understood): 



a) (°)(n|Di|n)(°) = Vi, 



^^\n\B^2\n) 



(0) 



b) ^°^(n|D.b)(°)^ "yi^r\ (°)(n|D.2by°) 



V2, 



E, 



(0) 



E. 



(0) 



£;. 



(0) 



E, 



(0) 



Wriy^j, 



(0) 



(n|^Ei|n)W = -(Vi^W), 



where F, 



F(x,), V, 



Dc, 



W(n|^Er|n)(°) = (V2i?(°)), 
Vj + igAj, and the transpose refers to the color 



matrices, we obtain at 0{l/m): 



n Q / ^ 



k^n 



(0) 



(fc|^Ei|n) 



(0) 



p(0) p(0) 



^(0,1) ^ 1 y 



kj^n 



(0) 



(A;|(7E^|n) 



'r|r.\(0) 



p(0) p(0) 



(20) 



By using translational invariance one can see that E^^'^^ and E^'^^ only depend on the 
relative distance r. Moreover, by using the symmetries of the static solutions, we can also 
see that E^'^^ = E^'^\ The expressions (|^) were first derived in Ref. 0. 
At 0(l/m^), we obtain 



c;(2,o) 



C 



(1)' 
^(0) 



.(1)2 



{n\[Bi.,gEi]\nr + --^Y: 



(0), 



n\gBi\k) 



(0) . (0) 



{k\gBi\n) 



(0) 



k^n 



E 



pIPi, 



(0) 



(n|^Ei|A;)W(o)(A;|c/Ei|n) 



(0) 



VlVj 



(o)(r2|^El|A;)(o)(o)(A;|^E>)W^ 



+2E 



(o)(n|^Ei|j)(o) {o)(^-|^E*^|A;)(o) (o)(A;|^E{|£)W (o)^£|^e>)W 



.(0) 



i5i°h3(^(0)_^(0))(^(0)_^(0). 



(o)(n|^Ei|A;)W (o)(A;|^Ei|j)W -W (j|^Ei|n)W 



+2 I vi E 

(o)(n|^Ei|A;)(o)(o)(fc|[Dr,^Ei]|n)(o) 



(i?l°)-i?f)3(E(°)-Ef) 



-VI- 



(i^;^^ 



(0) 



E^'y 



+3 VI- 



W(n|^El|A;)(o) W(A;|^Ei|n)W ■ {ViE^^^Y 



-2E 



(0) 



(n|(?Ei|j)(o) .(0) (j|(?Ei|A;)W W(A;|[Dr,5Ei]|n) 



(0) 



+6E 






+4. 



,W(n|[Dr,^Ei]|A:)W W(A:|^Ei|n)W ■ (Vi^j")) 



(i5i°) - 4°¥ 



(0) 



+ 



l W(n|[Dr,^Ei]|fc)WW(A:|[Dr,^Ei]|n) 

-4'^7a.c/rf'x^(°)(n|G»,(x)G;iJx)G:Jx)|n) 
(°)/ n|^Ei|A:)WW(A;|cri-^Bi|n) 



(0) 



.(1) 



+ 2 5 v.. 



(0) 



:l!hdM0) 



^o-i ■ (r X Vi 



(21) 



{E'n'-E'^'y J 4 r dr 

E(°'2)=e(2'0)(^Ei ^ -(7E^,(?Bi ^ -(?B^,o-i ^ 0-2, Vi ^ V2,Di ^ D,2,mi ^ ms), (22) 
and 

(o)(n|^El|A;)(o) (o)(A;|^En^)^°M 






v\vi 



(o)(n|^Ei|A;)(o)(o)(A;|^Enn)(°)^ 



-2E 



(0) _ r(0)^3 J 

(o)(n|^Ei|j)(o) W(j|^El|A;)(o) (o)(A;|^En^)^°^ (o)(£|^Enn)(°) 



+ |viE 



{e'^ - eF^e'^^ - Ef){E^^ - Ef\ 
(o)(n|(?Ei|A;)W W(A;|(7E^|j)(o) -^ (j|(7El^|n)(o)' 



(i?i°)-4°^)3(Ei°)-Ef)) 



- v^.E 



(0), 



(n|(?Ei|j)(").(")(j|(?Ei|A;)("HO)(^|^E^^i|^) 



{0)(0), 



^,i^i^\(o) 



,^n {Ei^^ - E^fiEi^^ - Ef) 

1 / , W(n|^El|A:)WW(fc|[D,2-,^El-]|n)W \ 

n^ (i?<?)-i?f)3 J 

1 / , W(0|[Dr,^Ei]|A:)WW(A:|(^Ef|n)W \ 
2^' (i5i°)-4°^)3 J 

_3 / . (o)(n|^E*i|A;)(o) (o)(A;|^E^|n)(o) ■ (V2^i°))\ 

"n ' (i?f - 4°V / 

_3( , (Vi^f ) -W (n|(7Ei|A;)W W(A;|(7Ef |n)W 

"2 r ^ (Ei°) - eF' . 



+ E 
-E 



(0) 



(^I^EiIj) 



(0) .(0) 



(j|^Ei|fc)(o)(o)(fc|[D,2-,^Ei^]|n) 



(0) 



{E'i^ ~ EfnE'i^ - Ef] 



(0) 



(n|[Dr,<7Ei]|fc)WW(fc|<?E^|j) 



i(0) .(0) 



(jI^e^I^) 



,Tlt,\(o) 






-3E 



(0) 



(n|^Ei|j)W -W (j|(7Ei|fc)W W(A;|^E^|n)W ■ (V^i^W) 



(i?i°) 



.(0)^ 



Er)\E)-. 



(0) 



E, 



{0)^ 



+3E 



ViEf ) -W (n|(?Ei|A;)(o) (o)(A;|(7E^|j)(o) -^ (j|(7E^|n)(o) 



(i^^ 



(0) 



Ef)4(E. 



(0) 



E. 



(0)^ 



+ - 



3W(n|[Dr,f7Ei]|A;)W(o)(A;|f7E^|n)(o) . (Vs^f ) 



.(0) 



^{0)^ 



^ro^ 



+ - 



2 (i?^ 

3(ViEf ) -W (n|(?Ei|A;)W W(j|[D,2-,^E^]|n)(°) 



(i^;^' 



(0) 



Efr 



^Vi^W) .(0) (n|(7Ei|A;)(o) W(A;|(7E^|n)(o) . (Vs^i^)) 



.(0) 



(0)n5 



(i5;r^ - E_ 

l W(n|[Dr,^Ei]|A:)WW(A:|[D,2-,^Ei^]|n)W 
r.iT|>,\(o)(o) 




(0) 



(nl^E^a |fc)WW(A:|cri-^Bi|n) 

(i?f - Efy 



(-En 



(0) 



E^ 



(0)^ 



J1)J2) 



E 



(0)/ 



n|o-i-^Bi|A;)WW(^|cr2.^B^|n) 



^lt.\(o) 



i?i°) - i?. 



(0) 



-(4.0-1 ■ 0-2 + d,J°)(n|TiVi ■ T2"^o-2|n)(°)) 5(3)(xi - X2) 



(23) 



The above equations (pOD-(p3D give the energies of the gluonic excitations between heavy 
quarks within an expansion in 1/m up to 0{l/m?). From these expressions, in the case of 
the ground state (ra = 0), we will derive, in section 0, the equivalent Wilson loop expressions. 

17| | in order to derive, from the QCD Hamil- 



A similar approach has been used in Ref. 



tonian in the Coulomb gauge, the spin-dependent part of the potential up to 0{l/m?). 
However, the behaviour at scales of 0{m) was not correctly incorporated there. If we take 
our NRQCD matching coefficients at tree level and neglect the tree-level annihilation con- 
tributions in the equal-mass case, we find agreement for the spin-dependent potentials (up 
to some transpose color matrices). Nevertheless, our general expression (18) differs from the 
one used in [0], which, in general, will not give the correct spin-independent potentials. 
This has to do, in our opinion, with the fact that in order to derive Eq. (^) one has to deal 
with operators rather than with numbers. 
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IV. PNRQCD 

In the previous section we have studied the static hmit of NRQCD and its corrections 
within a 1/m expansion. Let us now connect those results with pNRQCD. 

In the static limit, the gap between different states at fixed r will depend on the dimen- 
sionless parameter Aqcd^- In a general situation, there will be a set of states {n^s} such 
that E^^^(r) ~ mv'^ for the typical r of the actual physical system. We denote these states as 
ultrasoft. The aim of pNRQCD is to describe the behaviour of the ultrasoft states. There- 
fore, all the physical degrees of freedom with energies larger than mv'^ will be integrated 
out from NRQCD in order to obtain pNRQCD. It is in this context that one may work 
order by order in 1/m (in particular for the kinetic energy), and the calculation of the previ- 
ous section becomes the matching calculation between NRQCD and pNRQCD and provides 
a rigorous connection with the adiabatic approximation (this approximation is implicit in 
all the attempts at deriving the non-perturbative potentials from QCD we are aware of). 
Whereas this can be justified within a perturbative framework, in the non-perturbative case, 
we cannot, in general, guarantee the validity of the 1/m expansion and one may think of 
examples where certain degrees of freedom cannot be integrated out in the 1/m expansion 
p9[|). We believe that this possibility, which, to our knowledge, has never been men- 



see 



tioned before, except in Ref. [Q, deserves further study. Note that this does not have to do 
with the consideration of ultrasoft effects, which, unlike in earlier approaches, can be readily 
incorporated within our formalism. 

In the perturbative situation Aqcd'" -C 1, which has been studied in detail in [Q, {n^^s} 
corresponds to a heavy- quark-ant iquark state, in either a singlet or an octet configuration, 
plus gluons and light fermions, all of them with energies of 0{mv'^). In a non-perturbative 
situation, which we will generically denote by Aqcd^ ~ 1, it is not so clear what {n^s} 
is. One can think of different possibilities. Each of them will give, in principle, different 
predictions and, therefore, it should be possible to experimentally discriminate among them. 
In particular, one could consider the situation where, because of a mass gap in QCD, the 
energy splitting between the ground state and the first gluonic excitation is larger than 
mv"^, and, because of chiral symmetry breaking of QCD, Goldstone bosons (pions/kaons) 
appear. Hence, in this situation, {n^s} would be the ultrasoft excitations about the static 
ground state (i.e. the solutions of the corresponding Schrodinger equation), which will be 
named the singlet, plus the Goldstone bosons. If one switches off the light fermions (pure 
gluodynamics), only the singlet survives and pNRQCD reduces to a pure two-particle NR 
quantum-mechanical system, usually referred as a pure potential model. 

In this paper, we will study the pure singlet sector, with no reference to further ultrasoft 
degrees of freedom. In this situation, pNRQCD only describes the ultrasoft excitations about 
the static ground state of NRQCD. In terms of static NRQCD eigenstates, this means that 
only |0; xi,X2)*-°'' is kept as an explicit degree of freedom whereas |n; xi, X2)*-°^ with n 7^ 
are integrated outQ. This provides the only dynamical degree of freedom of the theory. 



^In fact, we are only integrating out states with energies larger than mv'^ and all the states 
with n 7^ will be understood in this way throughout the paper. Since, in practice, we are 
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It is described by means of a bilinear colour singlet field, S'(xi,X2,t), which has the same 
quantum numbers and transformation properties under symmetries as the static ground state 
of NRQCD in the one-quark-one-antiquark sector. In the above situation, the Lagrangian 
of pNRQCD reads 

^CpNRQCD = SUido - /is(xi,X2,Pi,P2)J5', (24) 

where hg is the Hamiltonian of the singlet (actually hg is only a function of r, pi, p2, which is 
analytic in the two last operators but typically contains non-analyticities in r), pi = — iVxi 
and p2 = — «Vx2- It has the following expansion up to order l/rn^: 



.2 „2 



^(1,0) ^(0,1) y(2,0) ^(0,2) yil,!) 



/i,(xi,X2,pi,p2) = ^ + ^ + \/(o) + ^— + ^— + ^-^ + ^-^ + ^— . (25) 
Zrrii 2m2 mi m2 rnf m2 mim2 

The integration of higher excitations is trivial using the basis |n; xi, X2) since, in this case, 
they are decoupled from |0; xi, X2). Then, the matching of NRQCD to pNRQCD consists in 
renaming things in a way such that pNRQCD reproduces the matrix elements of NRQCD 
for the ground state, and, in particular, the energy. This fixes the matching condition 

Eo(xi, X2, Pl, P2) = /is(xi, X2, pi, P2). (26) 

Although our main concern in this paper is to provide a well-controlled derivation of the 
potential for the heavy quarkonium, we would like to say a few words about the expressions 
En (n 7^ 0) we have found in the previous section. In the static limit, the different E^^ 
{n 7^ 0) are identified with the static potentials to be used in a Schrodinger equation to 
obtain the spectra of the bound systems composed of a heavy quark and an antiquark (plus 
glueballs) different from the heavy quarkonium such as, for instance, heavy hybrids. This 
assignment is argued within the adiabatic approximation and corresponds to what is actually 
done in lattice simulations [^. In this respect, since we have given a systematic method to 
obtain the corrections to the energy within a 1/m expansion, the energies En correspond to 
the quantum-mechanical Hamiltonians of the different bound systems made by a heavy quark 
and an antiquark (up to glueballs) and the 1/m and 1/m^ terms should be understood as the 
relativistic corrections to the static potentials. It is still an open problem if this procedure is 
the sensible thing to do for heavy hybrids, if (and whichever) other possibilities may occur, 
and if these potentials, like the heavy quarkonium potential, may eventually be written in 
terms of Wilson loops. We will not deal with these problems here, which, however, deserve 
further investigations. We refer to for related discussions. 



integrating over all the states, if we are in the situation where some states, different from the 
singlet, are ultrasoft, these have to be subtracted later on. This is analogous to what happens in the 
perturbative situation, where the subtraction is done order by order in the multipole expansion. In 
this situation our calculation should be understood as the leading term in the multipole expansion. 
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V. HEAVY QUARKONIUM POTENTIAL AND WILSON LOOPS 

In this section we express the heavy-quarkonium potential in terms of Wilson-loop op- 
erators. These kinds of expressions are quite convenient for lattice simulations or for QCD- 
vacuum- model studies (see, for instance, lHll^). We shall use the following definitions. The 



angular brackets (...) will stand for the average value over the Yang-Mills action, Wn for 
the rectangular static Wilson loop of dimensions r x T^: 



Wn = 'P exp I -ig i dz^'AJz 

I JrxTw 



and ((...)) = {. . .Wn)/ {Wn)] P is the path-ordering operator. Moreover, we define the 
connected Wilson loop with Oi(ti), 02(^2), ..., On(t„) operator insertions for Ti^/2 > ti > 
t2>--->tn> -Tw/2 by: 

((Oi(ti)02(t2)))e = ((Oi(ti)02(t2))) - ((Oi(t0))((O2(t2))), (27) 

((C'l(tl)02(t2)03(t3)))c = ((Oi(ti)02(t2)03(t3))) 

-((Oi(t0))((O2(t2)O3(t3)))c - {{0^it^)0,it,)))^{{Osih))) - {{0,ih))) {{0,(1,))) {{0,(1,))) , (28) 

((Oi(ti)02(t2)03(t3)04(t4)))c = {{Ol{h)02{t2)Os{h)0^{U))) 
-((Oi(ti)))((02(t2)03(t3)04(t4)))c-((Oi(ti)02(t2)))e((03(t3)04(t4)))c 
-((Oi(t0O2(t2)O3(t3)))c((O4(t4)))-((Oi(ti)))((O2(t2)))((O3(t3)O4(t4))). 
-((Oi(t0))((O2(t2)O3(t3))).((O4(t4)))-((Oi(t0O2(t2)))c((O3(t3)))((O4(t4))) 

-{{0,ih))){{02ih))){{0,it,))){{0,{U))), (29) 

We also define in a short-hand notation 

lim = lim lim , (30) 

where Tw is the time length of the Wilson loop and T the time length appearing in the time 
integrals. By performing first the Tw — *■ C)0, the averages ((. . .)) become independent of Tw 
and thus invariant under global time translations. 

By using the matching condition ( pG]) and the quantum- mechanical expressions (|20D, it 
has already been proved in |^ that the quarkonium singlet static potential and the Oll/m) 
potential can be expressed in terms of Wilson loops with field strength insertions in it as 

l-W(r)= hm Un{Wa), (31) 

I — *oo _/ 

V^''"\r) = -- hm f citt((^Ei(t)-(7Ei(0))),. (32) 

Owing to invariance under charge conjugation plus mi <-^ 1712 transformation we have 

\/(i.o)(r) = \/(0'i)(r). 
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The way to prove the equivalence of Eq. ( |32[ ) and Eq. ( pO[ ) has been discussed in Ref. [|], 
where more details can be found. Here we only mention that this equivalence proof as well as 
the following ones can be done straightforwardly by inserting complete sets of intermediate 
states in the Wilson loop operators and by explicitly computing the time integrals. 

Let us now consider the terms of 0(l/m^). It is convenient to split them in a spin- 
dependent and a spin-independent part. For the V^'^'^^ and V^^'"^^ potentials we define 

^(2.0) = kS^O) ^ ^^(2,0) ^ ^(0,2) ^ y(0,2) ^ ^J0,2) _ ^33^ 

The spin-independent terms can be written as 

J- r ..('2.nw x-i v/t2 r 



V^r' = -. {P?, V^''\r)} + ^^^^^Lj + V^^^\r), (34) 



and 



r 



1 T/(°'2)('^ 



VS'' = ^ [vl V$^\r)} + '-^^lA + Vr'\r), (35) 

where Li = r x pi and L2 = r x p2. Note that neither Li nor L2 corresponds to the 
orbital angular momentum of the particle and antiparticle. By using invariance under charge 
conjugation plus mi ^-> m2 transformation, we obtain 

Vy\r) = V^^'^\r), V^^'°)(r) = Vi^''\r), V}^^'\r) = V}'^^\r; m^ ^ mi). (36) 

The spin-dependent part of V"*^^'"-* is of the type 

^S°^ = vif)(r)Li-Si. (37) 

Analogously, for the V^^'"^"* potential we can write 

^i^'^ = -Vif (r)L2-S2. (38) 

From invariance under charge conjugation plus mi 4-^ 1712 transformation, we obtain 

Vi,f\r) = vi^^\r,m2^mi). 



By using Eqs. ( P6|) and (|2l|) we get, in terms of Wilson loop operators: 



vS''\r) = ^PF hm r dtf{{gE\it)gEim)^, (39) 



vS'"\r) = - (5^^ - 3FF) hm f dtt'{{gE^i{t)gEm))c. (40) 
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(1)' 
K(^'°)(r) = -^ lim {{[B„ gE,]{t)))c (41) 

• (1)2 ,, 

_^Ci^ K^ r ^.//.T3 m . .T3^m\\ , ^rT72T/(2.0)^ 



lim r dt{{gB,{t) ■ gBM))c + ^(V^V^^ 



4 T 
i rT rti pt2 



-'- lim / dt, f dt2 I ' rfts (^2 - t3)'((^Ei(ti) ■ f7Ei(t2)(7Ei(t3) ■ ^Ei(0)))e 
+]- fv; lim Trfti r rft2 (ti - t2)2((^El(ti)(7Ei(t2) ■ g^m))^ 
-I (V:y(°)) lim rdh r dh (ti - t2)'((^Ei(ti)(7Ei(t2) ■ ^Ei(O))), 

Z ^ ^ T^ooJq Jq 

--lim / rfti rdt2(ti-t2)'(([Di.,(7Ei](ti)(7Ei(t2)-(7Ei(0))), 

ZT^ooJo Jo 

+ ^ lim rrftt2^([Di.,(7Ei](t)[Di.,(?Ei](0)))e 
-'-(vilun^J^^dtt'{{gE\{t)[B,.,gE,]mc^ 
-llim rrftt3(([Di.,^Ei](t)^Ei(0))),(V^:yW) 

4 T^ooJo 

+\l^^l^^^[dtt'{{gE\it)gEm)),iViV^'^)^ 

--^lim rrftt^((^El(t)^Ei(0))),(V:yW)(Vfy(°)) 
-d?'fatc /ci^x lim ^((G;,(a;)Gj„(x)G:Jx))) 

(note that, although, formally the first and last terms depend on the time where the operator 
insertion is made, this is not so after doing the Tw -^ oo limitQ), 



vSf'ij) = -^ 



(1) rT ^(1) 



IV ■ lim / dtt{{gB^{t) X gE^)) + %^ ' (V,^^^)). (42) 



For the V^^^'^^ potential we define 

l^(M)=yjU)^^a,i)_ (43) 

The spin-independent part can be written as 

V^Y' = -I {Pi ■ P2, Vy\r)} - ^^^{U ■ L2 + L2 ■ L,) + V}''^\r), (44) 



2y(o) could also be written in a similar way: 

y(0) = - f d?^ lim (((n^n'' + B'^B") (x))). 

2 7 Tw^oo 
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while the spin-dependent part contains the following operators: 

^S'^ = v}^-^(r)U ■ S, - F£,\)(r)L2 ■ S^ + V^l''\r)S, ■ S, + V^l;'\r)S,,{r) , (45) 

where 812(1") = 3f • cri f • cr2 — cri • cr2. Because of the invariance under charge conjugation 
plus mi <-^ 1712 transformation, we have 

By using Eqs. ( p6|) and (|23|) we get, in terms of Wilson loop operators: 



Vy\r) = tr^r^ \im / dt f {{gE\{t) gEHO))) , 



(46) 



Z i— >ooJo 



(47) 



r2T/(l,l)^ 



rT rti /■t2 



(48) 



-I lim / dti r dt2 f ' dt3 (t2 - hf{{gEi{ti) ■ gEi{t2)gE2{h) ■ gE2m)c 

T^ooJo Jo Jo 

+ 1 fvt lim Trfti /*' dt2{h - t2y{{gE\{h)gE2{t2) ■ ^£2(0))),) 

2 y T^ooJo Jo J 

+ 1 fv; lim r^ti r dt2{h - t2y{{gEl{ti)gEi{t2) ■ gEi{0))))j 

2 \ T~*ooJo Jo J 

_1 (v;yW) hm Trfti r t/t2(ti - t2)'((^Ei(ti)^E2(t2) ■ ^E2(0))), 

2 ^ ' T~*oo Jo Jo 

_1 (v;y(°)) hm Tdti r rft2(tl - t2)'((^E^(tl)^Ei(t2) ■ ^Ei(0)))e 

--lim / C?ti /'rft2(tl-t2)'(([Di.,(7Ei](ti)^E2(t2)-(7E2(0)))e 

2 r^oo Jo Jo 
+ ^lim / rfti rrft2(ti-t2)'(([D2.,(7E2](ti)(7Ei(t2)-(7Ei(0)))e 

2 T^ooJo Jo 

-^ hm rc?tt2(([Di.,^Ei](t)[D2.,^E2](0))), 

4 T^oo Jo 

+ ^ (^V;^lim |^'^rftt2|^^^E.(^)P^_^^E2](0))),- ((^E^2W[Di-,^Ei]^ 

-^hm rrftt3{(([Di.,^Ei](t)^E^,(0)))e-(([D2.,^E2](t)^Ei(0))),}(V^,l^(°)) 

4 T^oo Jo '- J 

+ ^ [v;^lim /^"^rftt^ {((^El(t)^EKO)))e + ((^E^^W^EKO))),} (ViV^'^)) 



-- hm / dtt^(^Ei(t)<?EKO)))e(V;y(°))(V^,\/W) 

b T^oo Jo 

+ (4s + 4. ^hm {{T^T^))) <5(3)(xi - X2) 
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(here and in the following formulas the two colour matrices in {{T1T2)) are inserted in the 
Wilson loop at the same time: —Tw/2 < t < T]y/2; the t dependence disappears in the 
Tw — >■ cxD limit), 



V&iir) 






ir- lim / dtt {{gBi{t) x gEiiO))) 



(49) 



l/<,"'(r) 



9 (1) (2) T 

^^zlim / dt {{gB.it) . gB,m) 
-4(4. + 4. lim {{T^T^))) 5(=^)(xi - X2) 

T^ — s-oo 



(50) 



K 



(1,1)/ 



J1)J2) 



if*F lim 



dt 



((^Bl(t)^B^2(0)))-?^((^Bi(t)-^B2(0))) 



(51) 



We now compare our results with previous ones. For the spin-dependent potentials we 
find agreement with the Eichten-Feinberg results P] (once the NRQCD matching coefficients 
have been taken into account) except for the l/mim2 spin-orbit potential V^j^j • Since the 



Eichten-Feinberg results have been checked by, at least, three independent groups |rU, 12,131 
we perform a more detailed comparison in Appendix p. We show that our expression in terms 
of Wilson loops and theirs give different results in terms of intermediate states and, more 
important, we show that they give different perturbative results at leading order in as- Ours 
coincides with the well-known tree-level calculation, whereas the Eichten-Feinberg expression 
gives 1/2 the expected result. Moreover, our perturbative result fulfils the Gromes relation 
p!0| . The fact that the same mistake has been done by several groups can only be explained 



by a systematic error. We believe that their systematic error has to do with the common 
assumption in the literature that one may neglect, in general, the dependence of the Wilson 
loops on the gluonic strings, or on any other gluonic operator, at t = ±Tvi//2. An analysis of 
the calculation done by Eichten and Feinberg in supports this belief. Finally, we would like 
to mention that several different expressions for the spin-dependent potentials, in particular 
the correct one, can be found in the literature dealing with the lattice evaluation of them 
|30|-p2t^|l4| . All these refer to the work of Eichten and Feinberg [§ for the derivation. We 
believe that our result makes mandatory a clarification of all previous lattice evaluations of 
the spin-dependent potentials. 

The spin-independent potentials have only been computed before by Barchielli, Bram- 
billa, Montaldi and Prosperi in [jl2| (the analysis done in |Tl|, which appears to be incon- 
clusive, has never been published). We agree (once the NRQCD matching coefficients have 
been taken into account) with their results for the momentum-dependent terms, but not 
for the momentum-independent terms, where we find new contributions. Moreover, since 
the potential we get here is complete up to order 1/m^, it is not affected by the ordering 
ambiguity, which affects the derivation in [0. In this context, we would like to mention 
that our result may be of particular relevance for the study of the properties of the QCD 
vacuum in the presence of heavy sources. So far the lattice data for the spin-dependent and 
spin-independent potentials are consistent with a flux-tube picture, whereas it is only for the 
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spin-dependent terms that the so-called "scalar confinement" is consistent with lattice data 
|l|j33[| (however the lattice data are still not conclusive). It will be interesting to see how 
these pictures compare with the new momentum- and spin-independent potentials, once lat- 
tice data will be available for them. We note that some of them are not simply expressed by 
two field insertions on a static Wilson loop, such as the spin- or the momentum-dependent 
terms. In particular, an extended object coming from the Yang-Mills sector is required 
(similar extended objects would also show up by taking into account operators with light 
quarks) . 

A. Gauss law and further identities 

The above results may be simplified and rewritten in several ways. For instance, by using 
the quantum-mechanical identities a)-c) given in section |ITT|, we obtain 

hm (([Di,(7Ei](t))),= -fv^V(°)+22lim [^ dt {{gE,{t) ■ gE,m)j , (52) 

changing the expression of the Darwin term (that now looks similar to the analogous ex- 



pression given in Ref. |T2[)- In fact, by using the quantum-mechanical identities a)-c) of 
section |T|, we could systematically transform [D,5fE] in terms of normal derivatives acting 
on matrix elements or on static energies. 

Another possibility, which turns out to be more powerful, is the use of the Gauss law 
(P^). It allows us to write all the terms of the type [D, gfE] in terms of ^'•^•'(xi — X2) times 
some color matrices (up to some terms proportional to ^'^'^•'(O) that vanish in dimensional 
regularization). More information can be obtained by using the behaviour of the Wilson 
loops (or of the states) at short distances for the terms proportional to the deltas (assuming 
they are regular enough). It follows that all the original terms with [D, gE] disappear except 
the Darwin term. Moreover, we have (C/ = [N^ — l)/{2Nc)) 

hm ((T,%'^))5(=^)(xi - X2) = C/5(3)(xi - X2). 

Iw — *oo 

Therefore, some potentials get simplified into the following expressions 



yj2,0)(^) ^ ^C/^sCz, ^(3)(^^_^^) 



(53) 



hm / dt{{gB,it) ■ gB,{0)))c + ^{V^l 



'-^^ hm / ' " ^ 

4 T~*oo Jo 

- lim / dti / dt2 / 

2 T^oo Jo Jo Jo 

1 / . fT j-ti 



-'- lim / rfti / ' dt2 f ' dts (t2 - hfdgEiih) ■ gEi{t2)gEi{ts) ■ gEi{0))), 

ZT^ooJo Jo Jo 

+i fv; hm ^dti r dh (ti - t2)^{gE\{h)gE,it2) ■ ^Ei(0)))e') 

Z \ T^ooJo Jo J 

_1 (v;v(°)) hm ^dh r dt2 (ti - t2)^{gE\ih)gE^it2) ■ ^Ei(O))), 

Z ^ ' T^ooJo Jo 

+lhl^^^[dtt'{{gE\{t)gE{m),iViV^'^)\ 



--^lim r dtt\{gE\it)gEim)^iVlV^'^)iViV^'^) 
-d^'fatc fd'^ lim 9{{G;Ax)G'^^{x)GtM)) , 



V}'^'\r) = -\{ylV^^^) (54) 

-I lim / dtiT dh f ' dts {h - t3)'((^Ei(ti) • (7Ei(t2)^E2(t3) ■ ^E2(0)))e 

+ i fv; lim r dU r dhih - h)''{{gE\{t,)gE2{t2) ■ ^£2(0))),^ 

Z \ T^ooJo Jo J 

+ i fv; lim Trfti r rft2(ti - t2)'((^E^(ti)^Ei(t2) ■ 9E,{0)))c] 
_1 (vt\/(°)) lim Tdti r dt2(ti - t2)3(((7Ei(ii)^E2(t2) • gE^iO))), 

Z ^ ' T^oo Jo Jo 

_1 (Vt^w) lim Tdti r dt^ih - hf{{gEl{t,)gE,{h) ■ gE,{0))), 

2 V / T^ooJo Jo 

+ \ivl\imj\tt'{{{gE\{t)gEi{0)%+{{gEi{t)^^^^^ 

~ lim r dtt\{gE},{t)gEm))c{^lV^'^){^iV^'^) + (4. + d,,Cf) ^^^^xi - X2), 






lim / *{(jBi(()-jB2(0)))-4(<i,, + 4,C^)«<''>(x, -X2). (55) 



Similar considerations also apply to the results in terms of states of section 111. 



VI. POWER COUNTING 

The standard power counting of NRQCD (organized in powers oiv and as) used to assess 



the relative importance of the different matrix elements, as discussed, for instance, in [^ 
can only be proved in the perturbative regime. Even in this regime, owing to the different 
dynamical scales still involved, the matrix elements of NRQCD do not have a unique power 
counting in v. In the non-perturbative regime the problem of the power counting of NRQCD 
is still open. In principle, it is possible that a different power counting may be appropriate in 
this situation and this would influence, for instance, the studies of the charmonium system or 
of higher bottomonium states^. We believe that our result, through the connection between 
NRQCD and the quantum-mechanical picture, will eventually help to better understand 



^ A different, non-standard, power counting of the matrix elements of NRQCD may explain the 
apparent difficulties that NRQCD is facing to explain the polarization of prompt J/f/' data, and to 



accurately determine the different matrix elements (see [35|) 
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the hierarchy of the different matrix elements in NRQCD, as well as to get a much deeper 
understanding of the underlying dynamics. This is due to the fact that, by going to a NR 
quantum-mechanical formulation, we have made the dynamics of the heavy quarks explicit 
transfering the problem of the power counting of NRQCD into the problem of obtaining the 
power counting of the different potentials in pNRQCD. These may be expressed in terms of 
Wilson loops where only gluons and light quarks appear as dynamical entities and for which 
there are or there will be direct lattice measurements. Moreover, it is in this formulation that 
statements such as the virial theorem have a more rigorous, gauge-independent meaning. 

Here, we only say a few words about the expected behaviour of the potentials using 
arguments of naturalness on the scale mv, i.e. assuming that the potentials scale with mv. 
We first consider V^*^''). In principle, V^*^''^ counts as mv, but, by definition, the kinetic energy 
counts as mv"^. Therefore, the virial theorem constrains V^*^''^ also to count as mv'^. The 
extra 0{v) suppression has to come on dynamical grounds. In the perturbative case, it 
originates from the factor ag ^ v in the potential. In the non-perturbative case little can 
be said and other mechanisms must be responsible. Using naturalness, V^^'^"^ /m scales like 
mv^. Therefore, it could in principle be as large as V^^\ This makes a lattice calculation 
or a model evaluation of this potential urgent. Perturbatively, owing to the factor a^, it 
is of 0{mv'^). For what concerns the 1/m^ potentials, the naturalness argument suggests 
that they are of 0{mv^). However, also here several constraints apply. Terms involving 
Vy*^°'' ~ m^v^ are suppressed by an extra factor v, due to the virial theorem. The Gromes 



relation |T0|,|36 



+ Vj^r^ - V^ = 0, (56) 



^dV^'^ ^.(2,0) ,,(1,1) 
L2S1 

.(2,0) ..(1,1) 



2r dr ^^ ^^^ 



suppresses by an extra factor v the combination V^s ~ ^^2^1 • Similar constraints also exist 
for the spin-independent potentials [12|. Perturbatively the l/m? potentials count at most 



as 0{mv'^), because of the extra as suppression. Finally, it is important to consider that 
some of the potentials are 0(Q;s)-suppressed because of the matching coefficients inherited 
from NRQCD. This is, for instance, the case of the terms coming from the 1/m? corrections 
to the purely gluonic sector of the NRQCD Lagrangian or of the terms coming from the 
4-fermion sector. 

Terms involving two field-strength insertions in the static Wilson loop are known from 
lattice measurements ||1^ and have been studied in some QCD vacuum models |jT6|. For 



them a parameterization is possible and some supplementary information can be extracted. 
However, terms involving more than two field insertions in the static Wilson loop have not 
been studied so far, to our knowledge, by lattice simulations or within models. Consistency 
with the experimental data will further constrain any possible power-counting rule. In any 
case, a detailed study of the potentials using the above information (as well as new lattice 
or model-dependent results) should be performed in order to obtain the size (and thus 
the power-counting rules) of the different potentials for the charmonium and bottomonium 
systems. 
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VII. CONCLUSIONS AND OUTLOOK 

A new formalism with which to obtain the QCD potential at arbitrary orders in 1/m 
has been explained in detail. We have obtained expressions for the energies of the gluonic 
excitations between heavy quarks valid beyond perturbation theory at 0{l/m?). In par- 
ticular, for the heavy quarkonium, we have also obtained the complete spin-dependent and 
spin-independent potentials at 0{l/m?) for pure gluodynamics in terms of Wilson loops. 
For the spin-dependent piece our results correct the expressions given in [^,rU,r^|. For the 



spin-independent potentials, we agree with the momentum-dependent potentials obtained in 
I^ , but not for the momentum- independent terms, where new contributions are found. We 



have also briefly discussed the power counting in the non-perturbative regime. 

We conclude, commenting on two possible developments of the present work. First, it 
is worthwhile to explore the possibility of expressing the potentials associated with higher 
gluonic excitations in terms of Wilson loop operators as done here for the heavy quarkonium 
ground state. The corresponding quantum-mechanical expressions are given in Eqs. (pOl)- 
(^). Second, our results are complete at 0{l/rn?) in the case of pure gluodynamics. If we 



want to incorporate light fermions, the procedure to be followed is analogous and our results 
still remain valid (considering now matrix elements and Wilson loops with dynamical light 
fermions incorporated), except for new terms appearing in the energies at 0{l/m?) due to 
operators involving light fermions that appear in the NRQCD Lagrangian at 0{l/m?) |P4|| . 
They may be incorporated along the same lines as the terms discussed here and will be 
explicitly worked out elsewhere. 
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APPENDIX A: THE PARAMETERIZATION OF ||T^ 

In this appendix, for ease of comparison, we write the spin-independent potentials in the 
parameterization given in [|12|. They read 

V^f^ = \ jplpi, 5^^V,{r) + [^ - Ff) V,{r)\ + ^(V,2[F(°) + K]), (Al) 



Ki.i) 



'^ij 



V^r = 2 PlPi ^^'H(^) + Y - r^r' ^e(r) + Vf{r) 



(A2) 



Let us note that in |jT2| the l/mim2 potential contained only momentum-dependent pieces. 
Therefore, the momentum- independent potential, which we name V/, was missing. Our 



calculation also substantially modifies the result for Va given in |12|. The above potentials 
read 



yrf=^lim r dte{{g^^{t)-g¥.M))c 



(A3) 



'S^o 



SiJ 



Y - ^'^'j ^^ = 2 l^ool ^^^' ((^El(t)^Ei(O))), - -{{gE,{t) ■ ^Ei(O))), } , (A4) 



1 
8' 



V^[VW + K]) = K^^'°HO - \{^'rV^!''^) + \ (vyi\imJ^'^dtt'{{gE\{t)gEim)}j , 



(A5) 



H = -^ lim rdtf{{gE,{t)-gE,{0))), 



(A6) 



Y - f^fA K = -z\imj\tt' UgE\{t)gEi{0))), - ^-^{{gE^{t) ■ gE^iO)))} , (A7) 



r2T.{l,l)N , ^ 



y/(r) = V;^^'^^(r) + -(V,V^.^'^0 + o V^V^.li^ ciU^((^El(t)^E^2(0))) 



2 ^ " '■ ■ P^ ' ■ 2 V ' ' r^oo Jo 



(A8) 



APPENDIX B: COMPARISON WITH THE EICHTEN FEINBERG SPIN-ORBIT 

POTENTIAL 



(1) 



In order to compare our results with the Eichten-Feinberg ones properly, we set Cp 
p' = 1. Then, our Eq. (|i9| ) reads 



.(2) 
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VLls!ir) = -^r . lim f dttiigB,{t) x g^M)) (Bl) 



W(0|(?Bi|A;)W X W(A;|(7El^|0)(o) 



pert " 









+ 0(ai 



On the other hand, Eichten and Feinberg obtain (we actually use the expression in Minkowski 
space given in Ref. fl^): 

2^2^' Z^ 2^ g2 ,p(0) (O)wp(O) p(0)x 

_^ aoa^W(m|^Bi|fc)Wx W(fc|^E^|0)W 

2^2"^' 2^ 2^ nl fpiO) p{0)wp(0) p(0)x 

i ^ (o)(0|^Bi|fc)(o) X W(A;|^Ei^|0)(o) 






'^/'^s , ^/ 2 



where the a„(xi,X2) are defined by 

V'^(xi)0(xi,X2)x(x2)|vac) = ^a„(xi,X2)|n;xi,X2)^°\ 

n 

being 

0(y, x) = P exp l^tg J^ds (y - x) ■ A(x - s(x - y))| (B3) 

the end-point string used in the Wilson loop operators. Note that we have fixed Tw = T in 
Eq. ([B2| ), as corresponds to the procedure followed by Eichten and Feinberg. 

The above calculation makes manifest the disagreement of Eq. ( |BT| ) with Eq. ( P2| ) both 
at the perturbative level as well as in the representation in terms of intermediate states. A 
possible source of disagreement may be traced back in the original paper of Eichten and 
Feinberg |Q to their Eq. (4.9b), which seems to be incorrect. Finally, the reason of this last 
error seems to be the improper treatment of the Wilson loops in the large-time limit. 
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